Real space application of the mean-field description of spin glass dynamics 
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The out of equilibrium dynamics of finite dimensional spin glasses is considered from a point of view 
going beyond the standard 'mean-field theory' versus 'droplet picture' debate of the last decades. 
The main predictions of both theories concerning the spin glass dynamics are discussed. It is shown, 
in particular, that predictions originating from mean-field ideas concerning the violations of the 
fluctuation-dissipation theorem apply quantitatively, provided one properly takes into account the 
role of a spin glass coherence length, which plays a central role in the droplet picture. Dynamics in 
a uniform magnetic field is also briefly discussed. 
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Spin glasses have played a major and inspiring role 
in the rapidly developing field of the dynamics of glassy 
systems |Q. In particular, important developments have 
been achieved in the study of aging phenomena , which 
are encountered in several microscopic systems. In this 
context, two very different dynamic descriptions have 
emerged. The 'mean-field theory' consists in the exact 
solution of the dynamics of fully-connected (or equiva- 
lently infinite-dimensional) spin glass models The 
'droplet picture' is more phenomenological but directly 
addresses the problem of a real space (as opposed to the 
configurational space) description ||. 

It is a recurrent theme in the field to interpret numer- 
ical or experimental data as validating one description 
at the expense of the other Q. Technical tools of the 
dynamic mean-field theory have been argued to be nec- 
essary for a complete understanding of the aging of fi- 
nite dimensional spin glasses p]. In this Letter, we show 
that their range of potential validity is in fact wider than 
previously thought (and can apply even when no replica 
symmetry breaking is present), provided the theory is 
supplemented with the idea of a growing length scale, 
which we call the 'coherence length'. A growing equili- 
bration length, the 'domain size', is at the heart of the 
description of aging phenomena by the droplet picture. 
Although we cannot give a rigorous proof of their iden- 
tity, we argue that these two length scales have the same 
physical content. We propose a construction for the co- 
herence length and show numerically that it is propor- 
tional to the domain size. Last, we argue that the use 
of dynamical concepts coming from mean-field ideas does 
not necessarily imply the existence of a static replica sym- 
metry broken phase. 

The mean-field description of spin glass dynamics 
stems from the asymptotic solution of the dynamical 
equations for models which are statically solved by the 
Parisi replica symmetry breaking scheme Q. The be- 
havior of the system is encoded in the autocorrelation 
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function C(t,t w ) and the conjugated response function 
R(t,t w ). It is shown that the decay of the correlation 
involves a complex pattern of time scales organized in a 
hierarchical way. This 'dynamic ultramctricity' is a direct 
counterpart of the static Parisi solution for these mod- 
els j(3|. However, this feature is absent from all known ex- 
perimental and numerical data in three dimensions [QJ^] 
which show instead that the slow decay of the correlation 
(or the thermoremanent magnetization M(t, t w ) in exper- 
iments) is a one-time-scale process, C(t,t w ) ~ C(t/t w ) 
(or M(t, t w ) ~ M{t/t w )), for times t^t w . 

Non-trivial predictions are also made concerning the 
relation between R and C which satisfy at equilibrium 
the fluctuation-dissipation theorem (FDT), TR(t,t w ) = 
dt m C(t,t w ). A generalization of the FDT is obtained by 
introducing the function X(t,t w ) through M 



dt„ 



(1) 



with X(t,t w ) — 1 at equilibrium. In mean-field models, 
it can be shown that X(t,t w ) becomes at long times a 
single argument function, allowing the definition of the 
'fluctuation-dissipation ratio' (FDR) through M: 



x(q) = lim X(t, t w ) 

t,t„^oo C{t.t w )=q 



(2) 



Moreover, this purely out of equilibrium quantity is re- 
lated to the spin glass order parameter P{q), which 
measures the equilibrium distribution function of over- 
laps 



(q)= / dq'P(q'). 



(3) 



It has been further argued that Eqs. (Q)-(^) are true 
for finite dimensional glassy systems The existence 
of a FDR of the form (||) has been numerically in- 
vestigated in finite dimensional models [Qj through the 
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study of the physically accessible quantities C(t,t w ) and 
X(t,t w ) = // dt'R(t,t'). Eq. (|) is then graphically 
checked by representing the variations of x as a func- 
tion of C parameterized by the time difference t — t w , 
since Eq. (0) implies at large times 



x(t,t v 



1 

T 



dqx(q), 



(4) 



C(t,t w ) 



i.e. the obtained x(C) relation is independent of the 
time ||. At equilibrium, x = 1 and thus \ = (1 — C)/T. 
Numerically a x(C) parametric curve is obtained for a 
large £„, and compared to S(C, L) = J c dq dq'P(q', L) 
where P(q, L) is the Parisi function computed in a system 
of linear size L as large as possible. The coincidence of 
these curves was used to argue that both quantities were 
close to their limit, t Wl L — > oo, and to deduce informa- 
tions on the nature of the low-temperature phase [9 . 

A somewhat different picture has been put forward in 
Rcf. [|l0| and checked in the case of the 2D XY model: a 
real space view of the aging behavior as an equilibration 
process taking place on a growing 'coherence length' £(t w ) 
leads to a generalization of relations (^)-(Q) at finite times 
and sizes: 



T X (t,t w ) = S(C(t,t w ),£(t w )). 



(5) 



Eq. (^) states that the off- equilibrium properties of the in- 
finite aging system at finite-time t w are connected to the 
equilibrium properties of a system of finite-size £(t w ) jll] , 
independently of the t w , L — > oo limits. In other words, 
the system is quasi-equilibrated up to a length scale 
£(t w ). For the 2D XY model, the coherence length 
is thus proportional to the dynamic correlation length 

Inspired by dynamical mean-field theory, Equation (|5|) 
postulates the existence of a growing coherence length 
£{t w ) which represents the spatial extent to which the 
system is equilibrated at time t w . Interestingly, this is the 
basic assumption of the droplet picture that the spin glass 
dynamics is governed by large-scale excitations, whose 
size £ increases with time t w during the aging as £(t w ) ~ 
(lni^) 1 ^ @. Thus, we expect £{t w ) ~ £{t w ). 

Numerically, £,{t w ) can be extracted from a four-point 
correlation function |p|,p^|-^4[. Its relevance for the scal- 
ing properties of various physical quantities has been 
demonstrated in Refs. Experimentally, £(t w ) was 

extracted through the scaling behavior of the Zeeman 
energy Jll) . All these 3D investigations indicate a power 
law growth, £(t w ) ~ C (T) , with a(T) ~ 0.2T/T g . 

We now investigate numerically the validity of Eq. 
and then discuss its important consequences for the the- 
oretical description of spin glasses. Our aim is to demon- 
strate its wide range of applicability, independently of 
the existence of replica symmetry breaking, in the con- 
text of which it was first proposed p|§. For this purpose, 



we consider the low (but finite) temperature dynamics of 
the two-dimensional Edwards- Anderson model, for which 
the spin glass transition is at T c = p6| . At low tem- 
perature, the relaxation time becomes so large that the 
aging dynamics is very similar to the 3D case, including 
the power law growth £(t w ) ~ t^ 2T |l3|]. The model is 
defined by H = Jij s i s j-, where Sj (i = 1, • • • , N) 

are N = Lx L Ising spins located on the sites of a square 
lattice of linear size L. The sum runs over pairs 

of nearest neighbors. The Jij are random Gaussian vari- 
ables of mean and variance 1. 
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FIG. 1. Susceptibility - Correlation parametric curves ob- 
tained in the aging regime (circles) and from the static S(C, L) 
(full lines). The dashed line is the equilibrium FDT. (a) Con- 
stant waiting time, t w = 10 4 , T = 0.6, 0.5, 0.4, 0.3, 0.2 (from 
top to bottom), (b) Constant temperature, T = 0.4, and 
t w = 10 2 , 3.10 3 , 3.10 4 (from bottom to top). Equilibrium 
data: L = 8, 10, 18 (from bottom to top). 

To compute the autocorrelation function C(t,t w ) = 
N^ 1 J^i (si(t)si(t w )) ((•••) indicate an average over ini- 
tial conditions and • • • over realizations of the disor- 
der), and the susceptibility x(tytw), & very large system, 
L = 400, is quenched at the initial time t w = from 
a disordered state to T G [0.2,1.0]. The susceptibility 
x{t,t w ) is measured after applying a random magnetic 
field hi (taken from a Gaussian distribution of mean 
and variance ho 2 ) between times t w and t in each site. In 
the linear response regime (we used 0.02 < ho < 0.05), 
one gets x{t,t w ) = h ^ 2 N^ 1 Y /i {hiS i (t)). The equi- 
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librium P{q, L) is computed independently by equili- 
brating (using parallel tempering JT^]) samples of sizes 
L e [6, 24] and temperatures T <E [0.2, 1.2]. By definition, 
P(q, L) is the disorder-averaged histogram of the overlap 
q = N~ x sfs^ between two equilibrated copies (a, b) 
of the system. For the sizes and temperatures investi- 
gated, P(q, L) has its common 'non-trivial' structure [^), 
with a peak around a L-dependent value of the 'Edwards- 
Anderson parameter', and a tail extending towards q ~ 
values, although we know that lirnz^oo P{q, L) — 8(q) at 
all temperatures T > 0. 

We are now in position to compare the x(C) curves 
obtained in the aging situation, with S(C, L). Our results 
are summarized in Fig. [|. We show hrst in Fig. [j]-a the 
parametric curves for the same large waiting time, t w = 
10 4 , and different temperatures. The dynamic curves are 
qualitatively similar to the 3D case [to our knowledge, no 
such data are available in 2D]. At short times (1— C <C 1), 
the curves follow the equilibrium FDT relation, which 
they smoothly leave at longer times. 

The validity of Eq. @ is demonstrated in two steps. 
Fig. 1-a shows first that the dynamic curves follow, within 
our numerical precision, the curves of S(C, L) obtained 
by the double integration of P(q,L), for a given value 
of L. Note that the coincidence of these two indepen- 
dently computed functions on their whole support is a 
very strong requirement which unambiguously defines the 
coherence length, £{t w ) — L. The time evolution of £(t w ) 
is then followed for T = 0.4 in Fig. 1-b which shows 
that dynamic curves for increasing t w coincide with static 
curves for increasing sizes. For all temperatures investi- 
gated, we find a growth law for £(t w ) consistent with the 
2D growth laws reported in Ref. for £(t w ), showing 
that £(t w ) oc £(t w ). Physically, the coincidence between 
dynamic and static data means that at time t w , the sys- 
tem is locally equilibrated up to a coherence length £(t w ). 
The precise link between £(t w ) and £(t w ) is however a 
tricky point since the static P(q,£(t w )) is sensitive to the 
boundary conditions || . We use here (as is usual |J ) pe- 
riodic boundary conditions which leads to the estimate 
£{t w ) ~ 2£(t w ). The equality £{t w ) = £(t w ) would proba- 
bly be obtained computing P(q, L) in a box of size L in- 
side a much larger system, as proposed in Ref. [H . This 
is of course much more computationally demanding. 

At very low temperature (T < 0.3), we find that the 
agreement between statics and dynamics is not as good, 
with slightly different shapes for static and dynamic data. 
This is probably because even at very large times, the co- 
herence length is very small {£ ~ 1 — 3), so that the regime 
where Eq. (0) becomes valid is not reached. However, a 
similar trend has been observed in Ref. |l(J, attributed 
to the role of topological defects. No data are available 
yet in 3D or 4D spin glasses at very low temperature || , 
and this point should be checked. 

The dynamic data of Fig. |l|-a obtained at different tem- 
peratures can be collapsed using the scaling variables 



A = xT 1 ^ and B = (1 - C)^, with $ a free pa- 
rameter chosen to get the best collapse of the curves 
in Fig. |. For 1 - C < 1 (short times), the FDT im- 
plies B = A, while we obtain at large times the power 



law B 



A i-i/* ~ A o 



oG5 



i.e. 



2.3. This scal- 



ing behavior has been proposed in Ref. |j] as a dy- 
namic analog for the so-called Parisi- Toulouse approxi- 
mation used in a mean-field static context. The relation 
B ~ A 0A1 has been numerically found in 3D and 4D spin 
glasses, and compared to an expected mean-field behav- 
ior B ~ A - 5 ||. That the scaling works also in 2D 
explicitly shows that it is not necessarily connected to an 
underlying replica symmetry breaking, and weakens the 
evidence presented in other numerical works. 




(1 - C^- 
FIG. 2. Scaling behavior of the dynamic curves of Fig. ^j-a 
with $ ~ 2.3, guessed from the Parisi- Toulouse approxi- 
mation. The dashed lines are the relations B = A and 
B ~ A - 565 . 
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FIG. 3. Temperature dependence of the overlaps (q) and 
q-min, in a uniform magnetic field h — 0.4. 

We discuss now different important implications of the 
finite-time / finite-size connection described by Eq. (|^). 
From a pragmatic point of view, first, this relation im- 
plies that numerical studies of large aging systems or 
small equilibrated systems potentially contain the same 
informations. In our opinion, this fact has been largely 
underestimated in the spin glass literature, which often 
tries to overcome the difficulty of obtaining thermalized 
samples of large sizes by using large times in dynami- 
cal simulations. We exemplify this point by discussing 
the behavior of spin glasses in a uniform magnetic field. 
Various static tests of the existence of a spin glass phase 
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give inconclusive results (P(q, L) in a field differs from 
its mean-field shape, Binder cumulants do not cross in 
a field 0), while the dynamic behavior in a field has 
been claimed to clearly demonstrate the existence of a 
replica symmetry broken phase, in 3D and 4D pl|. This 
intriguing fact led us to perform in 2D the simulations 
of Refs. pj[ |. We have computed different values of the 
overlap, q m %n and (q) = J dq' P(q' ', L)q' which, accord- 
ing to the mean-field shape of P(q, L = oo) in a field, 
should be different jl9| ]; q m in is obtained dynamically 
as the infinite time extrapolation of the overlap between 
two independently aging copies of the system (T^]. Fig. || 
shows that, as in 3D and 4D, the two values differ at low 
enough T, Taken at face value, this result would lead to 
the erroneous conclusion that a replica symmetry break- 
ing transition occurs at T c (h — 0.4) ~ 0.5 [pofl . This 
shows that dynamic studies, using the same computer 
resources, are in fact probing the same length scales as 
equilibrium ones and therefore have to be taken with the 
same care. Much larger times (in dynamics) and sizes (in 
statics) should be used to show that there is no transition 
in a field in 2D. 

The validity of Eq. (JsJ) implies that the same excita- 
tions of length t(t w ) are governing the dynamics, i.e. the 
decay of the correlation as C(t,t w ) ~ C(£(t)/£(t w )), and 
are contributing to the low-q part of P(q,L), i.e. to the 
non-FDT part of the parametric curves in Fig. [j]. We em- 
phasize that Eq. (|5|) would not be satisfied at finite time 
during the coarsening in a pure ferromagnet (below T c ) 
where the domain walls governing the aging are absent in 
the static situation with usual boundary conditions. (See 
the related discussion of the role of vortices in Ref. |p0[.) 
Hence, we prefered the name 'coherence length' to 'do- 
main size', since these two names refer to qualitatively 
different physical situations. Our results give no indica- 
tion on the structure of the relevant excitations in spin 
glasses, but underlines the relevance (and the need) of 
finite-T equilibrium studies of large-scale excitations for 
the description of aging. It is not established, in particu- 
lar, that the description B in terms of compact droplets 
is correct 

From a fundamental point of view, our results give 
support to the picture of aging in spin glasses as the suc- 
cessive equilibration of excitations of increasing length 
scales, recently put forward to interpret temperature- 
cycling experiments f22|| . We speculate that the absence 
of £(t w ) in mean- field theory and hence of the result- 
ing multi- length-scale dynamics, is "compensated" by the 
dynamic ultrametricity |5| . This results in the inability of 
the theory of correctly predicting the simple t /t w -scalmg 
of dynamic functions. Note that the reported growth 
laws are inconsistent with the logarithmic law predicted 
by scaling arguments |3|,^2| , which would anyway lead to 
the unobserved logt/ log ^-scaling. A thorough investi- 
gation of this law seems necessary, and current experi- 
ments and simulations p3[ should clarify this point. 



Note finally that the extremely slow growth of the co- 
herence length implies that the length scales involved in 
the dynamics are relatively small, even at experimental 
times. The provocative idea that the thermodynamic 
limit might be of no practical importance directly follows. 
This also implies that Eq. (||) should apply in experi- 
ments, leading to results quantitatively similar to Fig. [I]. 
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